Damping of a nanomechanical oscillator strongly coupled to a quantum dot 
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We present theoretical and experimental results on the mechanical damping of an atomic force 
microscope cantilever strongly coupled to a self-assembled In As quantum dot. When the cantilever 
oscillation amplitude is large, its motion dominates the charge dynamics of the dot which in turn 
leads to nonlinear, amplitude-dependent damping of the cantilever. We observe highly asymmetric 
lineshapes of Coulomb blockade peaks in the damping that reflect the degeneracy of energy levels 
on the dot, in excellent agreement with our strong coupling theory. Furthermore, we predict that 
excited state spectroscopy is possible by studying the damping versus oscillation amplitude, in 
analogy to varying the amplitude of an ac gate voltage. 



Coupling a nanomechanical object to quantum elec- 
tronics provides a system that can be used to probe both 
the mechanics and the electronics with extreme sensitiv- 
ity. It has been predicted that the electronics may be 
used to measure the quantum nature of the mechanical 
object and the reverse — using the mechanics to mea- 
sure the quantum nature of mesoscopic electronics — was 
recently demonstrated with superconducting qubits 0. 
Electromechanical systems that have attracted consider- 
able attention recently include quantum shuttles H, and 
mechanics coupled to single electron transistors [4j, |5( or 
tunnel junctions In most systems studied both ex- 

perimentally and theoretically, the interaction between 
the electronic and mechanical components is weak. 

In this paper we study strong coupling effects, both 
theoretically and experimentally, in an electromechani- 
cal system consisting of a quantum dot capacitively cou- 
pled to an atomic force microscope (AFM) cantilever. 
Electrons tunneling on and off the dot effectively damp 
the cantilever, and this damping exhibits Coulomb block- 
ade peaks as a function of bias voltage similar to those 
well known in the dot conductance, even in the limit 
of weak coupling [1, 0, [l(J. It has long been predicted 
that level degeneracy on the dot leads to lineshape asym- 
metry of Coulomb blockade peaks in the conductance 
fl~0 ]. Recently, we observed corresponding temperature- 
dependent peak shifts in the damping at weak coupling 
[lfj , but the lineshape asymmetry was far too small to be 
measured before now. However, by driving the cantilever 
to large oscillation amplitudes we enter a regime of strong 
coupling where its motion strongly modifies the tunnel- 
ing rates on and off the dot, and leads to a dramatic 
enhancement of the lineshape asymmetry. This enhance- 
ment is much greater than expected from simply extrap- 
olating the weak coupling theory; it is a non-adiabatic 
effect that stems from the similarity of timescales for dy- 
namics of the cantilever and the dot. Furthermore, we 
predict that by measuring the damping versus bias volt- 
age and oscillation amplitude, strong coupling provides a 
means to perform excited state spectroscopy on the dot. 
Note that very different strong coupling effects unrelated 
to degeneracy were recently reported for a driven carbon 



nanotube coupled to an embedded dot [H, El • 

Our results show that AFM operated at large oscilla- 
tion amplitudes may be used to study degeneracy and 
level spacing, so-called shell structure, in quantum elec- 
tronic systems. The particular systems studied here, self- 
assembled quantum dots, are candidates for applications 
such as quantum information processing, and measur- 
ing their shell structure has attracted extensive research 
effort 14, 15, li|, 17 1. Our technique allows the level de- 
generacy to be read off from a single sweep of damping 
versus bias voltage, and offers the practical advantage 
that non-contact AFM is a way to address many dots 
one by one without the need for electrical contacts. 

Setup. — The mechanical oscillator is an AFM can- 
tilever with resonant frequency loq/2t: = 166 kHz, spring 
constant fco = 48 N/m, and intrinsic quality factor of 
typically Q ~ 2 x 10 5 . It is driven on resonance in self- 
oscillation mode at constant amplitude and mean tip- 
sample gap of 19 nm [1(3]. The cantilever is coated with 
a 10 nm Ti adhesion layer and a 20 nm Pt layer to ensure 
good electrical conductivity at low temperature. All data 
in this paper was collected at 5K. The sample is grown by 
chemical beam epitaxy, with the relevant features being 
uncapped InAs dots on top of a 20 nm InP tunnel bar- 
rier and a 10 nm InGaAs two-dimensional electron gas 
(2DEG) which acts as a back electrode. For full sample 
details see Ref. [l(|. The bias voltage Vb is applied to 
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FIG. 1: (a) Schematic of the setup. Electrons tunnel on and 
off the dot via a 2DEG back electrode, (b) Equivalent circuit 
diagram, where Ctip depends on x. 
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the 2DEG layer, with the cantilever electrically grounded 
(see Fig. [TJ . The potential drop between the 2DEG and 
the dot is cxVb, where a — C t i p /Cs is extracted from 
the experiment and Cs = Cy p + C2DEG is the total dot 
capacitance. The dot-cantilever coupling arises through 
Ctip, which depends on the tip position x. Electrons tun- 
nel between the 2DEG and the dot when Vb is sufficient 
to lift Coulomb blockade, while tunneling between the 
dot and cantilever is negligible due to the relatively large 
distance between them. The fluctuating charge on the 
dot results in both damping and a resonance frequency 
shift Auj of the cantilever; in the limit of weak coupling 
these are well described by linear response [la ]. Here 
we focus on the damping, which is provided in addition 
to the frequency shift by a phase-locked loop frequency 
detector and automatic gain controller [Toj . 

Model. — Provided the cantilever motion is small com- 
pared to the the tip-dot separation, we can assume that 
Ctip depends linearly on x and write the charging Hamil- 
tonian of the dot as 11911 



dependent rates Y + (T_) to add (remove) an electron, 
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where n is the number of electrons on the dot, j\f = 
— VeCtip/e is the dimensionless gate voltage (or con- 
trol charge), and Eq — e 2 /2Cs is the charging energy 
[2(ij ]. 7ic,o is the oscillator-independent part of Tic, and 
AH, osc , is a constant electrostatic modification of the 
oscillator potential. Interactions between the dot and 
oscillator are described by the final term with coupling 
strength A = - (2E c V B /e) (1 - a) dC tip /dx. We stress 
that the strong coupling effects discussed here occur de- 
spite Ctip remaining linear in x; they arise from the x- 
dependent tunneling rates discussed below. From Eq. (fl]) 
we see that the dot charge exerts a force An on the os- 
cillator; conversely, the oscillator position x changes the 
energy cost of adding or removing an electron on the dot. 

We focus on the bias range where or 1 extra electrons 
reside on the dot, with other charge states prohibited by 
Coulomb blockade; it is simple to generalize this to the 
case of n or n+1 electrons on the dot. The tunneling rates 
between the dot and back electrode are calculated from 
Fermi's golden rule accounting for the shell structure of 
the dot, i.e. the degeneracy of single particle levels. For 
a shell of degeneracy v occupied by n s h n electrons, there 
are ry+ = v — n s h n ways to add an electron, and once it 
has been added there are ry_ = n s hcii + 1 ways to remove 
it. The extra energy with 1 electron on the dot (i.e. 
the chemical potential difference between the dot and 
back electrode) is E(x) — 2E C (1/2 — A/") - Ax, which is 
modulated by x through the cou plin g. In the classical 
oscillator limit, Hujq <C ksT 2lj, [22j, this results in x- 
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where T is the tunneling rate to a single particle state and 
/ is the Fermi function. The asymmetry between adding 
and removing electrons is the root of the asymmetry in 
Coulomb blockade peaks [11|. 

We describe the coupled system using a master equa- 
tion for the charge on the dot combined with a Fokker- 
Planck equation for the phase space distribution of the 
oscillator [23|,[24|]. The central objects are the probabili- 
ties Pq(x, u) and P\{x, u) to find the oscillator at position 
x and velocity u with or 1 extra electrons on the dot; 
these satisfy master equations with x-dependent rates, 

d t P (x,u) = £ P + T-{x)P 1 - r+(x)P Q , (4) 
d t Pi{x,u) = APr+T+^Po -r_(z)Pi, (5) 

where C n = ujq (x — x n — T /fco) d u — ud x + jod u u de- 
scribes the evolution of a driven, damped harmonic os- 
cillator and x n = An/ kg is the equilibrium position with 
n electrons on the dot. The damping coefficient 70 is in- 
trinsic to the oscillator without coupling to the dot, and 
T is the external driving force. 

While it is straightforward to simulate the master 
equations directly, we gain further insight by focusing 
on the simpler dynamics of system averages. Following 
Ref. [Ill, we make the approximation that averages of 
products may be factorized into products of averages. 
While we lose correlations contained in Eqs. Q and ||5j), 
we find from comparison with full simulations that the 
asymmetric damping lineshape is still captured. Within 
this approximation we use Eqs. J3J and ([5]) to obtain 
coupled equations for the average quantities, 
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a , T + A(Pi) 



fco 



(x) - 70 (u) , 



a t (p 1 > = r+« a; »(Po>-r_((x»<p 1 >. 



(6) 
(7) 
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We seek a solution where the cantilever oscillates at 
constant amplitude a, replicating the experiment, such 
that (x(t)) — acos(woi). Ignoring the frequency shift 
due to tunneling (since Auj <C ujq) and assuming that 
the total damping is small (70 +71 <C too, justified 
self-consistently), we find that the effective, amplitude- 
dependent damping due to tunneling is given by [25} 



71 



nkoa 



2-jt/uq 



dt sin (u> t) (Px(t)) , 



(9) 



and obtain constant amplitude oscillations for drive 
T = fco/^o (7o + 7i) ( u )- Eq. © explicitly connects the 
damping to the time- varying dot charge, (n(t)) = (Pi(t)), 
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FIG. 2: (a) Calculated first Coulomb blockade damping peak 
versus bias voltage from simulation (dots) and semi-analytic 
theory (solid lines) for the oscillation amplitudes shown. The 
green dashed line is the linear response calculation. Gate volt- 
age is plotted in reverse for consistency with experiment, (b) 
Adiabatic approximation (dash-dotted) compared to semi- 
analytic theory (solid) and full simulation (dots) for a — 1 
nm. (c-e) Time dependence of average dot charge (solid) for 
a — 1 nm, at bias points marked in (a). Cantilever position is 
also shown (thin dashed) as a reference. We took 2Ec = 31 
meV, u/Y = 1 and A — 10 meV/nm. Other parameters are 
taken from the experiment (see Setup). 



specifically the part that is out of phase with the can- 
tilever position (x(t)). Note that damping arises even in 
the case of weak coupling, and is measurable using an os- 
cillator of sufficiently high quality factor [lfj , but Eq. © 
remains valid at strong coupling. It also reduces our cal- 
culation of 71 to solving Eq. ([5]) for (Pi(t)) numerically 
(inserting (x(t)} — acos(u>ot)), a much easier task than 
directly simulating Eqs. (g]) and ((SJ). 

Coulomb blockade peaks in the damping occur at 
charge degeneracy points, where the dot energy is equal 
with either or 1 electrons and charge fluctuations are 
maximal |lo| . Fig. Ufa) shows the first damping peak 
versus gate voltage for several oscillation amplitudes, cal- 
culated both using Eq. © and from direct simulation of 
Eqs. (@| and (0 following the approach of Ref. [24 1. 
We assume the level structure of a cylindrically symmet- 
ric dot, which includes a 2-fold degenerate s shell and 
a 4- fold degenerate p shell 26]. The simulated damp- 
ing is well described by linear response (green dashed) at 
weak coupling; note that even this peak is very slightly 
asymmetric as expected. As the oscillation amplitude is 
increased, the peak becomes broadened and highly asym- 
metric. The enhanced asymmetry at strong coupling is 
complctcy missed in an adiabatic approximation, where 
one assumes that the oscillator motion is much slower 
than tunneling (see Fig. [2](b)). This is not surprising: 
since ui ~ T, the cantilever can move significantly be- 
fore an electron tunnels on or off the dot. On the other 



hand, the damping calculated from our semi-analytic the- 
ory (see Eq. ([9])) agrees very well with the full simulation, 
so we use it to understand why the lineshape in the damp- 
ing is so highly asymmetric at strong coupling. 

Asymmetric lineshape. — The asymmetric lineshape of 
Coulomb blockade peaks is a result of the asymmetry 
between adding or removing electrons to or from a de- 
generate shell on the dot (cf. Eqs. ^ and ((3])). Con- 
sider the bias points c and e on either side of the peak 
in Fig. [2](a) , equal distances from its center such that 
the largest amplitude oscillator motion (broadest peak) 
is just large enough to swing N onto the charge degen- 
eracy point. One might guess that that a tunnel event 
near N = 1/2 is equally likely in both cases, but in fact 
tunneling is twice as likely to occur when starting from 
point e, where the dot is initially empty. This is because 
the rate to tunnel onto the dot near the charge degen- 
eracy point is r + ~ 2r/(0) = T (for the first peak in 
the 2-fold degenerate s shell), while the rate to tunnel 
off is only T_ - T(l - /(0)) = T/2. The asymmetry is 
apparent in the time dependence of (n(t)) at three bias 
points, shown in Fig. [DJc-e). Tunneling is more likely 
starting from point e, and this leads to increased damp- 
ing. Conversely, the situation is reversed in the second 
Coulomb blockade peak where there is only one way to 
add an electron to the half full s shell but two ways to re- 
move one once it is full. Thus, for large amplitudes when 
the asymmetry is visible, the lineshape provides a way 
to read off the shell degeneracy from a single Vb sweep: 
each peak is skewed away from the center of its shell. 

While a similar argument leads to a very slightly asym- 
metric lineshape at weak coupling, the asymmetry at 
strong coupling is much greater than we would expect 
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FIG. 3: (a) Experiment (solid) and theory (dashed) for the 
first Coulomb blockade damping peak at three oscillation am- 
plitudes. We converted eVu to E using a — 0.04 extracted at 
weak coupling. A single fit parameter value A = 7.8 meV/nm 
produced all three theory curves, (b) Experiment and theory 
for the second peak with A = 9.2 meV/nm. The peak is 
skewed in the opposite direction as the first, (c) Measured 
damping for the p shell; theory shown in (d) with A = 11 
meV/nm. Other parameters were taken from experiment as 
described in the text. 
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by extrapolation. This is because, for sufficiently strong 
coupling A, the change in gate voltage due to the os- 
cillator motion is greater than the thermal broadening 
from the Fermi distribution of electrons in the back elec- 
trode. In other words, the oscillator motion dominates 
over temperature, Aa > ksT. When this is satisfied, 
the harmonic distribution of the oscillator position P{x), 
peaked at the turning points of its motion, causes the 
most asymmetric tunneling rates, r±(±a) at the oscilla- 
tor extrema, to become especially important. The extra 
weighting of the most asymmetric rates leads to the dra- 
matically asymmetric lineshape at strong coupling. 

Finally in our discussion of the lineshape asymmetry, 
we point out the importance of relative timescalcs. For 
a slow oscillator, u> <C T, the adiabatic approximation is 
valid since the dot charge quickly equilibrates in response 
to the slow cantilever motion. In this case the damping is 
simply given by a weighted average of the linear response 
result taken over the oscillating gate voltage, and the 
lineshape asymmetry remains immeasurably small (see 
Fig-Hfb)). In the opposite limit, ui ^>T, the dot charge 
cannot respond to the rapid oscillator motion and damp- 
ing is suppressed. This can be seen from Eq. ([9]): for 
uj/T ^> 1, (Pi(i)) is roughly constant over one oscillator 
period and the damping becomes vanishingly small. The 
case that we have focussed on and measured is u> ~ T, 
where the interplay between the oscillation and tunnel- 
ing timescales leads to maximal and highly asymmetric 
damping. 

Measured Damping. — The experimentally measured 
cantilever damping is compared with theory in Fig. [3] 
In (a) and (b) we fit the first two Coulomb blockade 
peaks for three oscillation amplitudes (given in the leg- 
end) using Eq. ([9]). For each charge degeneracy point, we 
used a single fit parameter A to fit the peak at all three 
amplitudes, obtaining the values given in the caption. 
These are in good agreement with the values obtained 
from our weak coupling experiment on the same dot [Io| . 
We took the charging energy 2Eq =31 meV and lever 
arm a — 0.04 extracted at weak coupling. In princi- 
ple these parameters may also be found by fitting the 
strong coupling data directly, but we took advantage of 
our weak coupling results as a calibration and kept them 
fixed. Lastly, we point out that the cantilever damping is 
dominated by tunneling: here we find a peak value of up 
to 71/70 ~ 5, and at weak coupling (where the cantilever 
motion is small and does not move J\f off the charge de- 
generacy point) we measured as high as 71/70 ~ 20 [l(| • 

In Fig. Etc) and (d) we show the measured and theo- 
retical damping versus bias over the entire p shell. This is 
calculated by a straightforward extension of our deriva- 
tion of Eq. ([9]) to allow up to four electrons to occupy 
the 4-fold degenerate p shell. We find good qualitative 
agreement even in the crude approximation of constant 
charging energy, and using single values of A and T over 
the entire shell [27j]. We took 2E C = 20 meV for the 
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FIG. 4: Differential damping with respect to bias voltage, 
plotted versus bias (converted to E), and oscillation ampli- 
tude a. For sufficiently large amplitude, a third peak appears 
on the line E — Aa — A sp (white dashed line). Inset: cut along 
black dashed line. Parameters are the same as in Fig. [2] 



p shell (estimated from the peak spacing), and roughly 
aligned the peaks by adjusting the p shell level splitting 
phenomenologically. Once this was done, a single set 
of parameters was used to produce the damping spectra 
at all three amplitudes. Most importantly, in both the- 
ory and experiment the four peaks in the p-shell become 
five at large amplitudes, with peaks emerging directly be- 
tween the charge degeneracy points. This is completely 
consistent with our simple theory: at large amplitudes, 
the oscillator distribution P(x) is peaked at its extrema 
and contributions to tunneling are most important there. 
Thus, for amplitudes such that Aa — Eq, the tunneling 
is maximal when the bias voltage is at the midpoint be- 
tween two degeneracy points. 

Excited state spectroscopy. — Theoretically, our setup 
can be used to perform level spectroscopy on the dot, for 
example to measure the energy difference between the s 
and p shells, A sp . The oscillator is directly analogous 
to an ac gate voltage on the dot [28]. In the same way, 
when the oscillator motion is large enough to allow tran- 
sitions to multiple energy levels on the dot, the effective 
tunneling rate increases as does the damping. This is 
possible when the change in energy due to the oscillator 
is equal to the energy spacing, or Aa > A sp . At large 
amplitudes we expect a jump in 7 at the bias voltage 
where E = Aa — A sp . This leads to a peak in dj/dE 
that forms a line when plotted versus E and a, as seen 
in Fig. |U Measuring the slope and intercept of this line 
in experiment would directly provide A and A sp . 

Conclusions. — We have shown that the oscillation am- 
plitude is a useful new axis to exploit in using mechanical 
measurements to probe quantum electronic systems. We 
demonstrated our technique using self-assembled quan- 
tum dots; however, its implications extend to other quan- 
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turn electronic systems that can be placed on an insulat- 
ing surface with a back electrode. In particular, it should 
be possible to use AFM to study the level structure of sin- 
gle molecules. Due to the large spacing between energy 
levels, these may be studied at relatively high tempera- 
tures using further increased oscillation amplitudes. 

This work was supported by NSERC, FQRNT and CI- 
FAR. 
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